High-Throughput Spike Detection in Greenhouse Cultivated Grain Crops with Attention Mechanisms-Based Deep Learning Models.

Detection of spikes is the first important step toward image-based quantitative assessment of crop yield. However, spikes of grain plants occupy only a tiny fraction of the image area and often emerge in the middle of the mass of plant leaves that exhibit similar colors to spike regions. Consequently, accurate detection of grain spikes renders, in general, a non-trivial task even for advanced, state-of-the-art deep neural networks (DNNs). To improve pattern detection in spikes, we propose architectural changes to Faster-RCNN (FRCNN) by reducing feature extraction layers and introducing a global attention module. The performance of our extended FRCNN-A vs. conventional FRCNN was compared on images of different European wheat cultivars, including "difficult" bushy phenotypes from 2 different phenotyping facilities and optical setups. Our experimental results show that introduced architectural adaptations in FRCNN-A helped to improve spike detection accuracy in inner regions. The mean average precision (mAP) of FRCNN and FRCNN-A on inner spikes is 76.0% and 81.0%, respectively, while on the state-of-the-art detection DNNs, Swin Transformer mAP is 83.0%. As a lightweight network, FRCNN-A is faster than FRCNN and Swin Transformer on both baseline and augmented training datasets. On the FastGAN augmented dataset, FRCNN achieved a mAP of 84.24%, FRCNN-A attained a mAP of 85.0%, and the Swin Transformer achieved a mAP of 89.45%. The increase in mAP of DNNs on the augmented datasets is proportional to the amount of the IPK original and augmented images. Overall, this study indicates a superior performance of attention mechanisms-based deep learning models in detecting small and subtle features of grain spikes.

HiC-TE: a computational pipeline for Hi-C data analysis to study the role of repeat family interactions in the genome 3D organization.

MOTIVATION: The role of repetitive DNA in the 3D organization of the interphase nucleus is a subject of intensive study. In studies of 3D nucleus organization, mutual contacts of various loci can be identified by Hi-C sequencing. Typical analyses use binning of read pairs by location to reduce noise. We use binning by repeat families instead to make similar conclusions about repeat regions. RESULTS: To achieve this, we combined Hi-C data, reference genome data and tools for repeat analysis into a Nextflow pipeline identifying and quantifying the contacts of specific repeat families. As an output, our pipeline produces heatmaps showing contact frequency and circular diagrams visualizing repeat contact localization. Using our pipeline with tomato data, we revealed the preferential homotypic interactions of ribosomal DNA, centromeric satellites and some LTR retrotransposon families and, as expected, little contact between organellar and nuclear DNA elements. While the pipeline can be applied to any eukaryotic genome, results in plants provide better coverage, since the built-in TE-greedy-nester software only detects tandems and LTR retrotransposons. Other repeats can be fed via GFF3 files. This pipeline represents a novel and reproducible way to analyze the role of repetitive elements in the 3D organization of genomes. AVAILABILITY AND IMPLEMENTATION: https://gitlab.fi.muni.cz/lexa/hic-te/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Patients With Common Variable Immunodeficiency (CVID) Show Higher Gut Bacterial Diversity and Levels of Low-Abundance Genes Than the Healthy Housemates.

Common variable immunodeficiency (CVID) is a clinically and genetically heterogeneous disorder with inadequate antibody responses and low levels of immunoglobulins including IgA that is involved in the maintenance of the intestinal homeostasis. In this study, we analyzed the taxonomical and functional metagenome of the fecal microbiota and stool metabolome in a cohort of six CVID patients without gastroenterological symptomatology and their healthy housemates. The fecal microbiome of CVID patients contained higher numbers of bacterial species and altered abundance of thirty-four species. Hungatella hathewayi was frequent in CVID microbiome and absent in controls. Moreover, the CVID metagenome was enriched for low-abundance genes likely encoding nonessential functions, such as bacterial motility and metabolism of aromatic compounds. Metabolomics revealed dysregulation in several metabolic pathways, mostly associated with decreased levels of adenosine in CVID patients. Identified features have been consistently associated with CVID diagnosis across the patients with various immunological characteristics, length of treatment, and age. Taken together, this initial study revealed expansion of bacterial diversity in the host immunodeficient conditions and suggested several bacterial species and metabolites, which have potential to be diagnostic and/or prognostic CVID markers in the future.

Escherichia coli Strains Producing Selected Bacteriocins Inhibit Porcine Enterotoxigenic Escherichia coli (ETEC) under both In Vitro and In Vivo Conditions.

Enterotoxigenic Escherichia coli (ETEC) and Shiga toxin-producing E. coli (STEC) strains are the causative agents of severe foodborne diseases in both humans and animals. In this study, porcine pathogenic E. coli strains (n = 277) as well as porcine commensal strains (n = 188) were tested for their susceptibilities to 34 bacteriocin monoproducers to identify the most suitable bacteriocin types inhibiting porcine pathogens. Under in vitro conditions, the set of pathogenic E. coli strains was found to be significantly more susceptible to the majority of tested bacteriocins than commensal E. coli. Based on the production of bacteriocins with specific activity against pathogens, three potentially probiotic commensal E. coli strains of human origin were selected. These strains were found to be able to outcompete ETEC strains expressing F4 or F18 fimbriae in liquid culture and also decreased the severity and duration of diarrhea in piglets during experimental ETEC infection as well as pathogen numbers on the last day of in vivo experimentation. While the extents of the probiotic effect were different for each strain, the cocktail of all three strains showed the most pronounced beneficial effects, suggesting synergy between the tested E. coli strains. IMPORTANCE Increasing levels of antibiotic resistance among bacteria also increase the need for alternatives to conventional antibiotic treatment. Pathogenic Escherichia coli represents a major diarrheic infectious agent of piglets in their postweaning period; however, available measures to control these infections are limited. This study describes three novel E. coli strains producing antimicrobial compounds (bacteriocins) that actively inhibit a majority of toxigenic E. coli strains. The beneficial effect of three potentially probiotic E. coli strains was demonstrated under both in vitro and in vivo conditions. The novel probiotic candidates may be used as prophylaxis during piglets' postweaning period to overcome common infections caused by E. coli.

Letter to the Editor: Significant mutation enrichment in inverted repeat sites of new SARS-CoV-2 strains.

In a recently published paper, we have found that SARS-CoV-2 hot-spot mutations are significantly associated with inverted repeat loci and CG dinucleotides. However, fast-spreading strains with new mutations (so-called mink farm mutations, England mutations and Japan mutations) have been recently described. We used the new datasets to check the positioning of mutation sites in genomes of the new SARS-CoV-2 strains. Using an open-access Palindrome analyzer tool, we found mutations in these new strains to be significantly enriched in inverted repeat loci.

SARS-CoV-2 hot-spot mutations are significantly enriched within inverted repeats and CpG island loci.

SARS-CoV-2 is an intensively investigated virus from the order Nidovirales (Coronaviridae family) that causes COVID-19 disease in humans. Through enormous scientific effort, thousands of viral strains have been sequenced to date, thereby creating a strong background for deep bioinformatics studies of the SARS-CoV-2 genome. In this study, we inspected high-frequency mutations of SARS-CoV-2 and carried out systematic analyses of their overlay with inverted repeat (IR) loci and CpG islands. The main conclusion of our study is that SARS-CoV-2 hot-spot mutations are significantly enriched within both IRs and CpG island loci. This points to their role in genomic instability and may predict further mutational drive of the SARS-CoV-2 genome. Moreover, CpG islands are strongly enriched upstream from viral ORFs and thus could play important roles in transcription and the viral life cycle. We hypothesize that hypermethylation of these loci will decrease the transcription of viral ORFs and could therefore limit the progression of the disease.

TE-greedy-nester: structure-based detection of LTR retrotransposons and their nesting.

MOTIVATION: Transposable elements (TEs) in eukaryotes often get inserted into one another, forming sequences that become a complex mixture of full-length elements and their fragments. The reconstruction of full-length elements and the order in which they have been inserted is important for genome and transposon evolution studies. However, the accumulation of mutations and genome rearrangements over evolutionary time makes this process error-prone and decreases the efficiency of software aiming to recover all nested full-length TEs. RESULTS: We created software that uses a greedy recursive algorithm to mine increasingly fragmented copies of full-length LTR retrotransposons in assembled genomes and other sequence data. The software called TE-greedy-nester considers not only sequence similarity but also the structure of elements. This new tool was tested on a set of natural and synthetic sequences and its accuracy was compared to similar software. We found TE-greedy-nester to be superior in a number of parameters, namely computation time and full-length TE recovery in highly nested regions. AVAILABILITY AND IMPLEMENTATION: http://gitlab.fi.muni.cz/lexa/nested. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

What Can Long Terminal Repeats Tell Us About the Age of LTR Retrotransposons, Gene Conversion and Ectopic Recombination?

LTR retrotransposons constitute a significant part of plant genomes and their evolutionary dynamics play an important role in genome size changes. Current methods of LTR retrotransposon age estimation are based only on LTR (long terminal repeat) divergence. This has prompted us to analyze sequence similarity of LTRs in 25,144 LTR retrotransposons from fifteen plant species as well as formation of solo LTRs. We found that approximately one fourth of nested retrotransposons showed a higher LTR divergence than the pre-existing retrotransposons into which they had been inserted. Moreover, LTR similarity was correlated with LTR length. We propose that gene conversion can contribute to this phenomenon. Gene conversion prediction in LTRs showed potential converted regions in 25% of LTR pairs. Gene conversion was higher in species with smaller genomes while the proportion of solo LTRs did not change with genome size in analyzed species. The negative correlation between the extent of gene conversion and the abundance of solo LTRs suggests interference between gene conversion and ectopic recombination. Since such phenomena limit the traditional methods of LTR retrotransposon age estimation, we recommend an improved approach based on the exclusion of regions affected by gene conversion.

Correction to: Nested plant LTR retrotransposons target specific regions of other elements, while all LTR retrotransposons often target palindromes and nucleosome-occupied regions: in silico study.

[This corrects the article DOI: 10.1186/s13100-019-0186-z.].

pqsfinder web: G-quadruplex prediction using optimized pqsfinder algorithm.

MOTIVATION: G-quadruplex is a DNA or RNA form in which four guanine-rich regions are held together by base pairing between guanine nucleotides in coordination with potassium ions. G-quadruplexes are increasingly seen as a biologically important component of genomes. Their detection in vivo is problematic; however, sequencing and spectrometric techniques exist for their in vitro detection. We previously devised the pqsfinder algorithm for PQS identification, implemented it in C++ and published as an R/Bioconductor package. We looked for ways to optimize pqsfinder for faster and user-friendly sequence analysis. RESULTS: We identified two weak points where pqsfinder could be optimized. We modified the internals of the recursive algorithm to avoid matching and scoring many sub-optimal PQS conformations that are later discarded. To accommodate the needs of a broader range of users, we created a website for submission of sequence analysis jobs that does not require knowledge of R to use pqsfinder. AVAILABILITY AND IMPLEMENTATION: https://pqsfinder.fi.muni.cz, https://bioconductor.org/packages/pqsfinder. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Nested plant LTR retrotransposons target specific regions of other elements, while all LTR retrotransposons often target palindromes and nucleosome-occupied regions: in silico study.

BACKGROUND: Nesting is common in LTR retrotransposons, especially in large genomes containing a high number of elements. RESULTS: We analyzed 12 plant genomes and obtained 1491 pairs of nested and original (pre-existing) LTR retrotransposons. We systematically analyzed mutual nesting of individual LTR retrotransposons and found that certain families, more often belonging to the Ty3/gypsy than Ty1/copia superfamilies, showed a higher nesting frequency as well as a higher preference for older copies of the same family ("autoinsertions"). Nested LTR retrotransposons were preferentially located in the 3'UTR of other LTR retrotransposons, while coding and regulatory regions (LTRs) are not commonly targeted. Insertions displayed a weak preference for palindromes and were associated with a strong positional pattern of higher predicted nucleosome occupancy. Deviation from randomness in target site choice was also found in 13,983 non-nested plant LTR retrotransposons. CONCLUSIONS: We reveal that nesting of LTR retrotransposons is not random. Integration is correlated with sequence composition, secondary structure and the chromatin environment. Insertion into retrotransposon positions with a low negative impact on family fitness supports the concept of the genome being viewed as an ecosystem of various elements.

Bacterial but Not Fungal Gut Microbiota Alterations Are Associated With Common Variable Immunodeficiency (CVID) Phenotype.

Common Variable Immunodeficiency (CVID) is the most frequent symptomatic immune disorder characterized by reduced serum immunoglobulins. Patients often suffer from infectious and serious non-infectious complications which impact their life tremendously. The monogenic cause has been revealed in a minority of patients so far, indicating the role of multiple genes and environmental factors in CVID etiology. Using 16S and ITS rRNA amplicon sequencing, we analyzed the bacterial and fungal gut microbiota, respectively, in a group of 55 participants constituting of CVID patients and matched healthy controls including 16 case-control pairs living in the same household, to explore possible associations between gut microbiota composition and disease phenotype. We revealed less diverse and significantly altered bacterial but not fungal gut microbiota in CVID patients, which additionally appeared to be associated with a more severe disease phenotype. The factor of sharing the same household impacted both bacterial and fungal microbiome data significantly, although not as strongly as CVID diagnosis in bacterial assessment. Overall, our results suggest that gut bacterial microbiota is altered in CVID patients and may be one of the missing environmental drivers contributing to some of the symptoms and disease severity. Paired samples serving as controls will provide a better resolution between disease-related dysbiosis and other environmental confounders in future studies.

High-throughput analysis revealed mutations' diverging effects on SMN1 exon 7 splicing.

Splicing-affecting mutations can disrupt gene function by altering the transcript assembly. To ascertain splicing dysregulation principles, we modified a minigene assay for the parallel high-throughput evaluation of different mutations by next-generation sequencing. In our model system, all exonic and six intronic positions of the SMN1 gene's exon 7 were mutated to all possible nucleotide variants, which amounted to 180 unique single-nucleotide mutants and 470 double mutants. The mutations resulted in a wide range of splicing aberrations. Exonic splicing-affecting mutations resulted either in substantial exon skipping, supposedly driven by predicted exonic splicing silencer or cryptic donor splice site (5'ss) and de novo 5'ss strengthening and use. On the other hand, a single disruption of exonic splicing enhancer was not sufficient to cause major exon skipping, suggesting these elements can be substituted during exon recognition. While disrupting the acceptor splice site led only to exon skipping, some 5'ss mutations potentiated the use of three different cryptic 5'ss. Generally, single mutations supporting cryptic 5'ss use displayed better pre-mRNA/U1 snRNA duplex stability and increased splicing regulatory element strength across the original 5'ss. Analyzing double mutants supported the predominating splicing regulatory elements' effect, but U1 snRNA binding could contribute to the global balance of splicing isoforms. Based on these findings, we suggest that creating a new splicing enhancer across the mutated 5'ss can be one of the main factors driving cryptic 5'ss use.

The Impact of DNA Extraction Methods on Stool Bacterial and Fungal Microbiota Community Recovery.

Our understanding of human gut microbiota in health and disease depends on accurate and reproducible microbial data acquisition. The critical step in this process is to apply an appropriate methodology to extract microbial DNA, since biases introduced during the DNA extraction process may result in inaccurate microbial representation. In this study, we attempted to find a DNA extraction protocol which could be effectively used to analyze both the bacterial and fungal community. We evaluated the effect of five DNA extraction methods (QIAamp DNA Stool Mini Kit, PureLinkTM Microbiome DNA Purification Kit, ZR Fecal DNA MiniPrepTM Kit, NucleoSpin® DNA Stool Kit, and IHMS protocol Q) on bacterial and fungal gut microbiome recovery using (i) a defined system of germ-free mice feces spiked with bacterial or fungal strains, and (ii) non-spiked human feces. In our experimental setup, we confirmed that the examined methods significantly differed in efficiency and quality, which affected the identified stool microbiome composition. In addition, our results indicated that fungal DNA extraction might be prone to be affected by reagent/kit contamination, and thus an appropriate blank control should be included in mycobiome research. Overall, standardized IHMS protocol Q, recommended by the International Human Microbiome Consortium, performed the best when considering all the parameters analyzed, and thus could be applied not only in bacterial, but also in fungal microbiome research.

Quadruplex DNA in long terminal repeats in maize LTR retrotransposons inhibits the expression of a reporter gene in yeast.

BACKGROUND: Many studies have shown that guanine-rich DNA sequences form quadruplex structures (G4) in vitro but there is scarce evidence of guanine quadruplexes in vivo. The majority of potential quadruplex-forming sequences (PQS) are located in transposable elements (TEs), especially close to promoters within long terminal repeats of plant LTR retrotransposons. RESULTS: In order to test the potential effect of G4s on retrotransposon expression, we cloned the long terminal repeats of selected maize LTR retrotransposons upstream of the lacZ reporter gene and measured its transcription and translation in yeast. We found that G4s had an inhibitory effect on translation in vivo since "mutants" (where guanines were replaced by adenines in PQS) showed higher expression levels than wild-types. In parallel, we confirmed by circular dichroism measurements that the selected sequences can indeed adopt G4 conformation in vitro. Analysis of RNA-Seq of polyA RNA in maize seedlings grown in the presence of a G4-stabilizing ligand (NMM) showed both inhibitory as well as stimulatory effects on the transcription of LTR retrotransposons. CONCLUSIONS: Our results demonstrate that quadruplex DNA located within long terminal repeats of LTR retrotransposons can be formed in vivo and that it plays a regulatory role in the LTR retrotransposon life-cycle, thus also affecting genome dynamics.

pqsfinder: an exhaustive and imperfection-tolerant search tool for potential quadruplex-forming sequences in R.

MOTIVATION: G-quadruplexes (G4s) are one of the non-B DNA structures easily observed in vitro and assumed to form in vivo. The latest experiments with G4-specific antibodies and G4-unwinding helicase mutants confirm this conjecture. These four-stranded structures have also been shown to influence a range of molecular processes in cells. As G4s are intensively studied, it is often desirable to screen DNA sequences and pinpoint the precise locations where they might form. RESULTS: We describe and have tested a newly developed Bioconductor package for identifying potential quadruplex-forming sequences (PQS). The package is easy-to-use, flexible and customizable. It allows for sequence searches that accommodate possible divergences from the optimal G4 base composition. A novel aspect of our research was the creation and training (parametrization) of an advanced scoring model which resulted in increased precision compared to similar tools. We demonstrate that the algorithm behind the searches has a 96% accuracy on 392 currently known and experimentally observed G4 structures. We also carried out searches against the recent G4-seq data to verify how well we can identify the structures detected by that technology. The correlation with pqsfinder predictions was 0.622, higher than the correlation 0.491 obtained with the second best G4Hunter. AVAILABILITY AND IMPLEMENTATION: http://bioconductor.org/packages/pqsfinder/ This paper is based on pqsfinder-1.4.1. CONTACT: lexa@fi.muni.cz. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

p53 Specifically Binds Triplex DNA In Vitro and in Cells.

Triplex DNA is implicated in a wide range of biological activities, including regulation of gene expression and genomic instability leading to cancer. The tumor suppressor p53 is a central regulator of cell fate in response to different type of insults. Sequence and structure specific modes of DNA recognition are core attributes of the p53 protein. The focus of this work is the structure-specific binding of p53 to DNA containing triplex-forming sequences in vitro and in cells and the effect on p53-driven transcription. This is the first DNA binding study of full-length p53 and its deletion variants to both intermolecular and intramolecular T.A.T triplexes. We demonstrate that the interaction of p53 with intermolecular T.A.T triplex is comparable to the recognition of CTG-hairpin non-B DNA structure. Using deletion mutants we determined the C-terminal DNA binding domain of p53 to be crucial for triplex recognition. Furthermore, strong p53 recognition of intramolecular T.A.T triplexes (H-DNA), stabilized by negative superhelicity in plasmid DNA, was detected by competition and immunoprecipitation experiments, and visualized by AFM. Moreover, chromatin immunoprecipitation revealed p53 binding T.A.T forming sequence in vivo. Enhanced reporter transactivation by p53 on insertion of triplex forming sequence into plasmid with p53 consensus sequence was observed by luciferase reporter assays. In-silico scan of human regulatory regions for the simultaneous presence of both consensus sequence and T.A.T motifs identified a set of candidate p53 target genes and p53-dependent activation of several of them (ABCG5, ENOX1, INSR, MCC, NFAT5) was confirmed by RT-qPCR. Our results show that T.A.T triplex comprises a new class of p53 binding sites targeted by p53 in a DNA structure-dependent mode in vitro and in cells. The contribution of p53 DNA structure-dependent binding to the regulation of transcription is discussed.

Effects of Isoflavone-Enriched Feed on the Rumen Microbiota in Dairy Cows.

In this study, we compared the effects of two diets containing different isoflavone concentrations on the isoflavone transfer from feed into milk and on the rumen microbiota in lactating dairy cows. The on-farm experiment was conducted on twelve lactating Czech Fleckvieh x Holstein cows divided into two groups, each with similar mean milk yield. Twice daily, cows were individually fed a diet based on maize silage, meadow hay and supplemental mixture. Control group (CTRL) received the basal diet while the experimental group (EXP) received the basal diet supplemented with 40% soybean isoflavone extract. The average daily isoflavone intake in the EXP group (16 g/day) was twice as high as that in the CTRL group (8.4 g/day, P<0.001). Total isoflavone concentrations in milk from the CTRL and EXP groups were 96.89 and 276.07 µg/L, respectively (P<0.001). Equol concentrations in milk increased from 77.78 µg/L in the CTRL group to 186.30 µg/L in the EXP group (P<0.001). The V3-4 region of bacterial 16S rRNA genes was used for metagenomic analysis of the rumen microbiome. The experimental cows exhibited fewer OTUs at a distance level of 0.03 compared to control cows (P<0.05) and reduced microbial richness compared to control cows based on the calculated Inverse Simpson and Shannon indices. Non-metric multidimensional scaling analysis showed that the major contributor to separation between the experimental and control groups were changes in the representation of bacteria belonging to the phyla Bacteroidetes, Proteobacteria, Firmicutes, and Planctomycetes. Surprisingly, a statistically significant positive correlation was found only between isoflavones and the phyla Burkholderiales (r = 0.65, P<0.05) and unclassified Betaproteobacteria (r = 0.58, P<0.05). Previous mouse and human studies of isoflavone effects on the composition of gastrointestinal microbial populations generally report similar findings.

Hammock: a hidden Markov model-based peptide clustering algorithm to identify protein-interaction consensus motifs in large datasets.

MOTIVATION: Proteins often recognize their interaction partners on the basis of short linear motifs located in disordered regions on proteins' surface. Experimental techniques that study such motifs use short peptides to mimic the structural properties of interacting proteins. Continued development of these methods allows for large-scale screening, resulting in vast amounts of peptide sequences, potentially containing information on multiple protein-protein interactions. Processing of such datasets is a complex but essential task for large-scale studies investigating protein-protein interactions. RESULTS: The software tool presented in this article is able to rapidly identify multiple clusters of sequences carrying shared specificity motifs in massive datasets from various sources and generate multiple sequence alignments of identified clusters. The method was applied on a previously published smaller dataset containing distinct classes of ligands for SH3 domains, as well as on a new, an order of magnitude larger dataset containing epitopes for several monoclonal antibodies. The software successfully identified clusters of sequences mimicking epitopes of antibody targets, as well as secondary clusters revealing that the antibodies accept some deviations from original epitope sequences. Another test indicates that processing of even much larger datasets is computationally feasible. AVAILABILITY AND IMPLEMENTATION: Hammock is published under GNU GPL v. 3 license and is freely available as a standalone program (from http://www.recamo.cz/en/software/hammock-cluster-peptides/) or as a tool for the Galaxy toolbox (from https://toolshed.g2.bx.psu.edu/view/hammock/hammock). The source code can be downloaded from https://github.com/hammock-dev/hammock/releases. CONTACT: muller@mou.cz SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Transposable elements and G-quadruplexes.

A significant part of eukaryotic genomes is formed by transposable elements (TEs) containing not only genes but also regulatory sequences. Some of the regulatory sequences located within TEs can form secondary structures like hairpins or three-stranded (triplex DNA) and four-stranded (quadruplex DNA) conformations. This review focuses on recent evidence showing that G-quadruplex-forming sequences in particular are often present in specific parts of TEs in plants and humans. We discuss the potential role of these structures in the TE life cycle as well as the impact of G-quadruplexes on replication, transcription, translation, chromatin status, and recombination. The aim of this review is to emphasize that TEs may serve as vehicles for the genomic spread of G-quadruplexes. These non-canonical DNA structures and their conformational switches may constitute another regulatory system that, together with small and long non-coding RNA molecules and proteins, contribute to the complex cellular network resulting in the large diversity of eukaryotes.